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O ■ INTRODUCTION 

7— I ■ 

^ ■ Non-perturbative aspects are believed to play a crucial role in understanding the 

ON ■ formation of a mass gap in the spectrum of excitations in a non-Abelian gauge theory, 

but despite much progress a simple explanation is still lacking. Bere this problem will 
be addressed in a finite volume, where its size can be used as a control parameter, which 
is absent in infinite volumes. The starting point, as in solving a quantum mechanical 
i-Ch ■ problem, is finding a proper description of the classical potential and its minima. 

>• ■ Due to the non-Abelian nature of the theory not only the classical potential, but 

also the kinetic term in the Lagrangian or Bamiltonian, is more complicated. The latter 
is a manifestation of the non-trivial Riemannian geometry of the physical configuration 
space 1 , formed by the set of gauge orbits A/Q (A is the collection of connections, Q 
the group of local gauge transformations). Most frequently, coordinates of this orbit 
space are chosen by picking a representative gauge field on the orbit in a smooth and 
preferably unique way. It is by now well known that linear gauge conditions like the 
Landau or Coulomb gauge suffer from Gribov ambiguities 2 . The reason behind this 
is that topological obstructions prevent one from introducing affine coordinates 3 in a 
global way. In principle therefore, one can introduce different coordinate patches with 
transition functions to circumvent this problem 4 . One way to make this specific is to 
base the coordinate patches on the choice of a background gauge condition 5, 6 . One 
could envisage to associate to each coordinate patch ghost fields and extend BRST 
symmetry to include fields with non-trivial "Grassmannian sections", although such a 
formulation is still in its infancy. 

We will pursue, however, the issue of finding a fundamental domain for non-Abelian 
gauge theories 7 and its consequence for the glueball spectrum in intermediate volumes. 
The finite volume context allows us to make reliable statements on the non-perturbative 
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contributions, because asymptotic freedom guarantees that at small volumes the effec- 
tive coupling constant is small, such that high-momentum states can be treated per- 
turbatively. Only the lowest (typically zero or near-zero momentum) states will be 
affected by non-perturbative corrections. We emphasize that it is essential that gauge 
invariance is implemented properly at all stages. We will describe the results mainly 
in the context of a Hamiltonian picture 8 with wave functionals on configuration space. 
Although rather cumbersome from a perturbative point of view, where the covariant 
path integral approach of Feynman is vastly superior, it provides more intuition on 
how to deal with non-perturbative contributions to observables that do not vanish in 
perturbation theory. An essential feature of the non-perturbative behaviour is that the 
wave functional spreads out in configuration space to become sensitive to its non-trivial 
geometry. If wave functionals are localised within regions much smaller than the in- 
verse curvature of the field space, the curvature has no effect on the wave functionals. 
At the other extreme, if the configuration space has non-contractible circles, the wave 
functionals are drastically affected by the geometry, or topology, when their support 
extends over the entire circle. Instantons are of course the most important examples of 
this. Not only the vacuum energy is affected by these instantons, but also the low-lying 
glueball states and this is what we are after to describe accurately, albeit in sufficiently 
small volumes. The geometry of the finite volume, to be considered here, is the one 
of a three-torus 9, 5 and a three-sphere 10 ' 11 . These lecture notes are an updated and 
extended version of ref 12 . 



COMPLETE GAUGE FIXING 

An (almost) unique representative of the gauge orbit is found by minimising the 
L 2 norm of the vector potential along the gauge orbit 7 ' 13 

F A {g) = ||M|| 2 = - J^x tr ^-^9 + g-%g) 2 ) , (1) 

where the vector potential is taken anti-hermitian. Expanding around the minimum of 
eq. (0), writing g(x) = exp(X(x)), one easily finds: 

\\ 9 A\\ 2 = \\A\\ 2 + 2 [ tr(Xd i A i )+ [ tr (X^FP(A)X) 

JM JM 

+\f tr(X[[A i ,X],d i X}) + ^- [ ti([D l X,X}[d t X,X]) + 0(x 5 ).(2) 

6 JM Yl JM v ' 

Where FP(A) is the Faddeev-Popov operator (ad(A)X = [A, X]) 

FP(A) = -diDi(A) = -d'i-OMiA,). (3) 

At any local minimum the vector potential is therefore transverse, d{Ai = 0, and 
FP(A) is a positive operator. The set of all these vector potentials is by definition the 
Gribov region Q. Using the fact that FP(A) is linear in A, Q is seen to be a convex 
subspace of the set of transverse connections T. Its boundary dQ is called the Gribov 
horizon. At the Gribov horizon, the lowest eigenvalue of the Faddeev-Popov operator 
vanishes, and points on dfl are hence associated with coordinate singularities. Any 
point on dQ can be seen to have a finite distance to the origin of field space and in 
some cases even uniform bounds can be derived 14 ' 15 . 

The Gribov region is the set of local minima of the norm functional (3) and needs 
to be further restricted to the absolute minima to form a fundamental domain, which 
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will be denoted by A. The fundamental domain is clearly contained within the Gribov 
region. To show that also A is convex we note that 

rA\\ 2 -\\A\\ 2 = I tr (A 2 )- J tr ([g- l A i9 + g'^gf) 

= Jtr (g*FP f (A) g) = (g, FP f (A) g) , (4) 

where FPf(A) is the Faddeev-Popov operator generalised to the fundamental represen- 
tation. At the critical points A E T of the norm functional (recall T = {A E A\diAi = 
0}) FPf(A), like FP(A), is a hermitian operator. We can define A in terms of the 
absolute minima over g E Q of (g, FPj(A) g) 

A = {A E T\ min (g, FP f (A) g) = 0}. (5) 

Using that FPf(A) is linear in A and assuming that A^ and A^ are in A and therefore 
satisfy the equation min q (g,FP f (A) g) = 0, we find that A = sA^ + (1 - s)A^ 
satisfies the same identity for all s E [0, 1] (such that both s and (1 — s) are positive). 
The line connecting two points in A, therefore lies within A. 

If we would not specify anything further, as a convex space is contractible, the 
fundamental region could never reproduce the non-trivial topology of the configuration 
space. This means that A should have a boundary 16 . Indeed, as A is contained in 
Q, this means A is also bounded in each direction. Clearly A = is in the interior 
of A, which allows us to consider a ray extending out from the origin into a given 
direction, where it will have to cross the boundary of A and Q. For any point along 
this ray in A, the norm functional is at its absolute minimum as a function of the 
gauge orbit. However, for points in Q that are not also in A, the norm functional is 
necessarily at a relative minimum. The absolute minimum for this orbit is an element 
of A, but in general not along the ray. Continuity therefore tells us that at some point 
along the ray, this absolute minimum has to pass the local minimum. At the point 
they are exactly degenerate, there are two gauge equivalent vector potentials with the 
same norm, both at the absolute minimum. As in the interior the norm functional has 
a unique minimum, again by continuity, these two degenerate configurations have to 
both lie on the boundary of A. This is the generic situation. 

It is important to note that A = is a so-called reducible connection 17 which has 
a non-trivial stabiliser, i.e. a subgroup of the gauge group that leaves the connection 
invariant, in this case the set of constant gauge transformations. These reducible con- 
nections give rise to so-called orbifold singularities, that manifest themselves through 
curvature singularities. We know perfectly well how to deal with it by not fixing the 
constant gauge transformations. Indeed the norm functional is degenerate along the 
constant gauge transformations and the Coulomb gauge does not fix this gauge degree 
of freedom. We simply demand that the wave functional is in the singlet representa- 
tion under the constant gauge transformations. The degeneracy of the norm functional 
along the constant gauge transformations gives rise to a trivial part in the kernel of 
the Faddeev-Popov operator (also directly seen by inspection, since it vanishes when 
acting on constant Lie-algebra elements that generate the constant gauge transforma- 
tions). This trivial kernel is to be projected out, something that is possible in a finite 
volume. In this formalism there is a remnant gauge invariance G, which requires no 
gauge fixing since its volume is finite, but still needs to be divided out to get the proper 
identification 

A/G = A/g. (6) 
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Here A is assumed to include the non-trivial boundary identifications. It is these 
boundary identifications that restore the non-trivial topology of AjQ. 

If the degeneracy at the boundary is continuous along non-trivial directions one 
necessarily has at least one non-trivial zero eigenvalue for FP(A) and the Gribov hori- 
zon will touch the boundary of the fundamental domain at these so-called singular 
boundary points. We sketch the general situation in figure 1. By singular we mean 
here a coordinate singularity. In principle, by choosing a different gauge fixing in the 
neighbourhood of these points one could resolve the singularity. If singular bound- 
ary points would not exist, all that would have been required is to complement the 
Hamiltonian in the Coulomb gauge with the appropriate boundary conditions in field 
space. Since the boundary identifications are by gauge transformations the boundary 
condition on the wave functionals is simply that they are identical under the boundary 
identifications, possibly up to a phase in case the gauge transformation is homotopically 
non-trivial, as will be discussed further on in more detail. 

Unfortunately, the existence of non-contractible spheres 3 in the configuration space 
allows one to argue that singular boundary points are to be expected 16 . Consider the 
intersection of a rf-dimensional non-contractible sphere with A. The part in the interior 
of A is contractible and it can only become non-contractible through the boundary iden- 
tifications. The simplest way this can occur 16 is if for the (<i-l)-dimensional intersection 
with the boundary, all points are to be identified. It would imply degeneracy of the the 
norm functional on a (rf-l)-dimensional subspace, leading to at least d-1 zero-modes for 
the Faddeev-Popov operator. The intersection with the boundary of A can, however, 
exist of more than one connected component. In the case of two such components 
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(c) 

Figure 2. Sketch of a singular boundary point due to a bifurcation of the norm func- 
tional. It can be used to show that there are homotopically trivial gauge copies inside 
the Gribov horizon (h). 

one can make a non-contractible sphere by identifying points of the (d-l)-dimensional 
boundary intersection of the first connected component with that of the second and 
there is no necessity for a continuous degeneracy. Conversely, not all singular boundary 
points, even those associated with continuous degeneracies, need to be associated with 
non-contractible spheres. 

When a singular boundary point is not associated to a continuous degeneracy, the 
norm functional undergoes a bifurcation moving from inside to outside the fundamental 
(and Gribov) region. The absolute minimum turns into a saddle point and two local 
minima appear, as indicated in figure 2. These are necessarily gauge copies of each 
other. The gauge transformation is homotopically trivial as it reduces to the identity 
at the bifurcation point, evolving continuously from there on. For reducible connec- 
tions, that have a non-trivial stabiliser, this argument may be false 18 , but examples of 
bifurcations at irreducible connections were explicitly found for S 3 , see ref. 19 (app. A). 
We will come back to this. 

Also Gribov's original arguments for the existence of gauge copies 2 (showing that 
points just outside the horizon are gauge copies of points just inside) can be easily 
understood from the perspective of bifurcations in the norm functional. It describes 
the generic case where the zero-mode of the Faddeev-Popov operator arises because of 
the coalescence of a local minimum with a saddle point with only one unstable direction. 
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At the Gribov horizon the norm functional locally behaves in that case as X 3 , with X 
the relevant zero eigenfunction of the Faddeev-Popov operator. The situation sketched 
in figure 2 corresponds to the case where the leading behaviour is like X 4 . See ref. 16 for 
more details and a discussion of the Morse theory aspects that simplify the bifurcation 
analysis. As the Gribov region is associated with the local minima, and since the space 
of gauge transformations resembles that of a spin model, the analogy with spin glasses 
makes it unreasonable to expect that the Gribov region is free of further gauge copies. 
This will be illustrated by explicit examples. Unfortunately restrictions to a subset 
of the transverse gauge fields is a rather non-local procedure. This cannot be avoided 
since it reflects the non-trivial topology of field space. 

GAUGE FIELDS ON THE THREE-TORUS 

Homotopical non-trivial gauge transformations are in one to one correspondence 
with non-contractible loops in configuration space, which give rise to conserved quan- 
tum numbers. The quantum numbers are like the Bloch momenta in a periodic potential 
and label representations of the homotopy group of gauge transformations. On the fun- 
damental domain the non-contractible loops arise through identifications of boundary 
points (as will be demonstrated quite explicitly for the torus in the zero-momentum 
sector). Although slightly more hidden, the fundamental domain will therefore contain 
all the information relevant for the topological quantum numbers. Sufficiently accurate 
knowledge of the boundary identifications will allow for an efficient and natural projec- 
tion on the various superselection sectors (i.e. by choosing the appropriate "Bloch wave 
functionals"). All these features were at the heart of the finite volume analysis on the 
torus 5 and we see that they can in principle be naturally extended to the full theory, 
thereby including the desired 9 dependence. In the next section this will be discussed 
in the context of the three-sphere. In ref. 6 we proposed formulating the Hamiltonian 
theory on coordinate patches, with homotopically non-trivial gauge transformations 
as transition functions. Working with boundary conditions on the boundary of the 
fundamental domain is easily seen to be equivalent and conceptually much simpler to 
formulate. If there would be no singular boundary points this would have provided a 
Hamiltonian formulation where all topologically non-trivial information can be encoded 
in the boundary conditions. Still, for the low-lying states in a finite volume, both on the 
three-torus and the three-sphere, singular boundary points will not play an important 
role in intermediate volumes. 

Probably the most simple example to illustrate the relevance of the fundamental 
domain is provided by gauge fields on the torus in the abelian zero-momentum sector. 
For definiteness let us take G = SU{2) and = i^Jz (L is the size of the torus). 
These modes are dynamically motivated as they form the set of gauge fields on which the 
classical potential vanishes. It is called the vacuum valley (sometimes also referred to as 
toron valley) and one can attempt to perform a Born-Oppenheimer-like approximation 
for deriving an effective Hamiltonian in terms of these "slow" degrees of freedom. To 
find the Gribov horizon, one easily verifies that the part of the spectrum for FP(A) 
that depends on C, is given by \ 9 ^{C) = 2nn- (2nn±C), with n ^ an integer vector. 
The lowest eigenvalue therefore vanishes if Ck = ±2n. The Gribov region is therefore 
a cube with sides of length 4n, centred at the origin, specified by |C fc | < 27r for all k, 
see figure 3. 

The gauge transformation g^) = exp(7iixfcT3/L) maps Ck to Ck + 27r, leaving the 
other components of C untouched. As gt^) is anti-periodic it is homotopically non- 
trivial (they are 't Hooft's twisted gauge transformations 20 ). We thus see explicitly 
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Figure 3. A two dimensional slice of the vacuum valley along the (Cj, C2) plane. The 
fat square give the Gribov horizon, the grey square is the fundamental domain. The 
dots at the Gribov horizon are Gribov copies of the origin. 



that gauge copies occur inside Q, but furthermore the naive vacuum A = has (many) 
gauge copies under these shifts of 27r that lie on the Gribov horizon. It can actually 
be shown for the Coulomb gauge that for any three-manifold, any Gribov copy by a 
homotopically non-trivial gauge transformation of A = will have vanishing Faddeev- 
Popov determinant 16 . Taking the symmetry under homotopically non-trivial gauge 
transformations properly into account is crucial for describing the non-perturbative 
dynamics and one sees that the singularity of the Hamiltonian at Gribov copies of 
A = 0, where the wave functionals are in a sense maximal, could form a severe obstacle 
in obtaining reliable results. 

To find the boundary of the fundamental domain we note that the gauge copies C = 
(ir, C 2 , C 3 ) and C = (— ir, C 2 , C 3 ) have equal norm. The boundary of the fundamental 
domain, restricted to the vacuum valley formed by the abelian zero-momentum gauge 
fields, therefore occurs where \Ck\ = it, well inside the Gribov region, see figure 3. The 
boundary identifications are by the homotopically non-trivial gauge transformations 

The fundamental domain, described by |C&| < tt, with all boundary points regular, 
has the topology of a torus. To be more precise, as the remnant of the constant gauge 
transformations (the Weyl group) changes C to — C, the fundamental domain A/G 
restricted to the abelian constant modes is the orbifold T 3 /Z 2 . Generalisations to 
arbitrary gauge groups were considered in ref. 6 . (The fundamental domain turns out 
to coincide with the unit cell or "minimal" coordinate patch defined in ref. 6 ). 

Formulating the Hamiltonian on A, with the boundary identifications implied 
by the gauge transformations g^), avoids the singularities at the Gribov copies of 
A = 0. "Bloch momenta" associated to the 2ir shift, implemented by the non-trivial 
homotopy of g^k), label 't Hooft's electric flux quantum numbers 20 ^{Ck = — fl") = 
exp(7uefc)\l/(Cfc = 7r). Note that the phase factor is not arbitrary, but ±1. This is 
because g^ is homotopically trivial. In other words, the homotopy group of these anti- 
periodic gauge transformations is Z\. Considering a slice of A can obscure some of the 
topological features. A loop that winds around the slice twice is contractible in A as 
soon as it is allowed to leave the slice. Indeed including the lowest modes transverse to 
this slice will make the Z 2 nature of the relevant homotopy group evident 5 . It should be 
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mentioned that for the torus in the presence of fields in the fundamental representation 
(quarks), only periodic gauge transformations are allowed. In that case it is easily seen 
that the intersection of the fundamental domain with the constant abelian gauge fields 
is given by the domain | C fc | < 2n, whose boundary coincides with the Gribov horizon. 
It is interesting to note that points on dQ form an explicit example of a continuous 
degeneracy due to a non-contractible sphere 15 . 

In weak coupling Liischer 9 showed unambiguously that the wave functionals are 
localised around A = 0, that they are normalisable and that the spectrum is discrete. 
In this limit the spectrum is insensitive to the boundary identifications (giving rise to 
a degeneracy in the topological quantum numbers). This is manifested by a vanishing 
electric flux energy, defined by the difference in energy of a state with \e\ — 1 and the 
vacuum state with e = 0. Although there is no classical potential barrier to achieve this 
suppression, it comes about by a quantum induced barrier, in strength down by two 
powers of the coupling constant. This gives a suppression 21 with a factor exp(— S/g) 
instead of the usual factor of exp(— 8n 2 /g 2 ) for instantons 22 . Here S = 12.4637 is the 
action computed from the effective potential. At stronger coupling the wave functional 
spreads out over the vacuum valley and the boundary conditions drastically change the 
spectrum 5 . At this point the energy of electric flux suddenly switches on. 

Integrating out the non-zero momentum degrees of freedom, for which Bloch de- 
generate perturbation theory provides a rigorous framework 23, 9 , one finds an effective 
Hamiltonian. Near A = 0, due to the quartic nature of the potential energy {F"-) 2 for 
the zero-momentum modes (the derivatives vanish and the field strength is quadratic 
in the field), there is no separation in time scales between the abelian and non-abelian 
modes. Away from A = one could further reduce the dynamics to one along the vac- 
uum valley, but near the origin this would be a singular decomposition (the adiabatic 
approximation breaks down). However, as long as the coupling constant is not too 
large, the wave functional can be reduced to a wave function on the vacuum valley near 
dA where the boundary conditions can be implemented. These boundary conditions are 
formulated in a manner that preserves the invariance under constant gauge transfor- 
mation and the effective Hamiltonian is solved by Rayleigh-Ritz (providing also lower 
bounds from the second moment of the Hamiltonian). The influence of the boundary 
conditions on the low-lying glueball states is felt as soon as the volume is bigger than 
an inverse scalar glueball mass. We summarise below the ingredients that enter the 
calculations. 

The effective Hamiltonian is expressed in terms of the coordinates c", where % = 
{1, 2, 3} is the spatial index (c = 0) and a = {1, 2, 3} is the SU(2)-colour index. These 
coordinates are related to the zero- momentum gauge fields through A®(x) = cf/L. 
We note that the field strength is given by F^ = —e a M(\c d /L 2 and we introduce the 
gauge-invariant "radial" coordinate r\ = y/^2 a cfcf. The latter will play a crucial role in 
specifying the boundary conditions. For dimensional reasons the effective Hamiltonian 
is proportional to 1/L. It will furthermore depend on L through the renormalised 
coupling constant {g{L)) at the scale n — 1/L. To one-loop order one has (for small 
L) g(L) 2 = 12tt 2 /[—ll ln(A ms L)}. One expresses the masses and the size of the finite 
volume in dimensionless quantities, like mass-ratios and the parameter z = mL. In 
this way, the explicit dependence of g on L is irrelevant. This is also the preferred way 
of comparing results obtained within different regularisation schemes (i.e. dimensional 
and lattice regularisation). The effective Hamiltonian is now given by 
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+71 Y r i + 72 Y r t + 73 Y r i r 1 + 74 Y r i + ^Y r i r 1 + 7e n r i 

i i i>j i ij^j i 

+a 3 Y r l F !k + «4 Y r\F% + « 5 det 2 c. (7) 

ijk,a ij,a 

We have organised the terms according to the importance of their contributions, 
ignoring terms quartic in the momenta. The first line gives (when ignoring ai j2 ) the 
lowest order effective Hamiltonian, whose energy eigenvalues are 0(g 2 ^ 3 ), as can be 
seen by rescaling c with g 2 ^ 3 . Thus, in a perturbative expansion c = 0(g 2 ^ 3 ). The 
second line includes the vacuum-valley effective potential (i.e. the part that does not 
vanish on the set of abelian configurations). These two lines are sufficient to obtain the 
mass-ratios to an accuracy of better than 5%. The third line gives terms of 0(g 4 ) in the 
effective potential, that vanish along the vacuum-valley. The coefficients (to two-loop 
order for 7^) are 
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The choice of boundary conditions, associated to each of the irreducible represen- 
tations of the cubic group 0(3, 7L) and to the states that carry electric flux 20 , is best 
described by observing that the cubic group is the semidirect product of the group of 
coordinate permutations S3 and the group of coordinate reflections Z\. We denote the 
parity under the coordinate reflection c" — > — c" by Pi = ±1 (i fixed). The electric flux 
quantum number for the same direction will be denoted by qi = ±1. This is related 
to the more usual additive (mod 2) quantum number by qj = exp(i7re J ). Note that 
for SU(2) electric flux is invariant under coordinate reflections. If not all of the electric 
fluxes are identical, the cubic group is broken to S2 x Z 3 , where S , 2(~ Z 2 ) corresponds 
to interchanging the two directions with identical electric flux (unequal to the other 
electric flux). If all the electric fluxes are equal, the wave functions are irreducible rep- 
resentations of the cubic group. These are the four singlets A^ 2 y which are completely 
(anti-)symmetric with respect to S 3 and have each of the parities pi = ±1. Then there 
are two doublets E ± , also with each of the parities pi — ±1 and finally one has four 
triplets T^Ly Each of these triplet states can be decomposed into eigenstates of the 
coordinate reflections. Explicitly, for Tf^ we have one state that is (anti-)symmetric 
under interchanging the two- and three-directions, with p 2 = p 3 = —p\ = +1. The 
other two states are obtained through cyclic permutation of the coordinates. Thus, 
any eigenfunction of the effective Hamiltonian with specific electric flux quantum num- 
bers qi can be chosen to be an eigenstate of the parity operators pi. The boundary 
conditions of these eigenfunctions ^^(c) are simply given by 

A( r .^( c ))| rj _ = , if Pt q t = +1 (9) 

and one easily shows that with these boundary conditions the Hamiltonian is hermitian 
with respect to the innerproduct < ty' >= f ri<w d 9 c^*(c)^'(c). For negative parity 
states (TiiPi = —1) this description is, however, not accurate 24 as parity restricted to 



9 



m' / m 




z=mL 



Figure 4. Mass ratios m'/m as a function of L in units of the inverse scalar glueball 
mass m _1 , z = mL. The analytic results are given by the full (continuum) and dashed 
(4 3 lattice) curves. Represented are the the square root of the energy of electric flux 
per unit length (bottom), the E tensor mass (middle) and the T2 tensor mass (top). 
The latter two are (almost) degenerate below z = 0.95, which is also where the electric 
flux energy is exponentially suppressed. Where no error bars on the Monte Carlo data 
are visible they are smaller than the size of the data points. 

the vacuum valley is equivalent to a Weyl reflection (a remnant of the invariance under 
constant gauge transformations). 

After correcting for lattice artefacts 25 , the (semi-) analytic results agree extremely 
well with the best lattice data 26 (with statistical errors of 2% to 3%) up to a volume 
of about .75 fermi, or about five times the inverse scalar glueball mass. In figure 4 we 
present the comparison for a lattice of spatial size 4 3 . Monte Carlo data 26 are most 
accurate for this lattice size. For more detailed comparisons see ref. 25 . The analytic 
results below z = 0.95 are due to Liischer and Miinster 27 , which is where the spectrum 
is insensitive to the identifications at the boundary of A. 

Most conspicuously the tensor state in finite volumes is split in a doublet E, with a 
mass that is roughly 0.9 times the scalar A 1 mass and a triplet T 2 with a mass of roughly 
1.7 times the scalar mass. Note that the multiplicity weighted average is approximately 
1.4 times the scalar mass, agreeing well with what was found at large volumes from 
lattice data 26 . 

Apart from the corrections for the lattice artefacts, generalisation to SU(3) was 
established by Vohwinkel 28 , with qualitatively similar results. In large volumes the rota- 
tional symmetry should be restored, as is observed from lattice simulations. The prop- 
erties of the fundamental domain restricted to the zero-momentum modes for SU(3) 
can be read off from the results in ref. 6 . In this reference also the generalisation to 
arbitrary gauge groups is discussed. 

At large volumes extra degrees of freedom start to behave non-perturbatively. To 
demonstrate this, the minimal barrier height that separates two vacuum valleys that 
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are related by gauge transformations with non-trivial winding number 



^'^/m^ 1 ^- (10) 

was found to be 72.605/ Lg 2 , using the lattice approximation and carefully taking the 
continuum limit 29 . As long as the states under consideration have energies below this 
value, the transitions over this barrier can be neglected and the zero-momentum effec- 
tive Hamiltonian provides an accurate description. One can now easily find for which 
volume the energy of the level that determines the glueball mass (defined by the differ- 
ence with the groundstate energy) starts to be of the order of this barrier height. This 
turns out to be the case for L roughly 5 to 6 times the correlation length set by the 
scalar glueball mass. The situation is sketched in figure 5. 

We expect, as will be shown for the three-sphere, that the boundary of the fun- 
damental domain along the path in field space across the barrier (which corresponds 
to the instanton path if we parametrise this path by Euclidean time t), occurs at the 
saddle point (which we call a finite volume sphaleron) in between the two minima. The 
degrees of freedom along this tunnelling path go outside of the space of zero-momentum 
gauge fields and if the energy of a state flows over the barrier, its wave functional will 
no longer be exponentially suppressed below the barrier and will in particular be non- 
negligible at the boundary of the fundamental domain. Boundary identifications in this 
direction of field space now become dynamically important too. The relevant "Bloch 
momentum" is in this case obviously the 8 parameter, as wave functionals pick up a 
phase factor e %e under a gauge transformation with winding number one. For many 
of the intricacies in describing instantons on a torus we refer to ref. 30 ' 31, 32 . On the 
three-torus we have therefore achieved a self-contained picture of the low-lying glueball 
spectrum in intermediate volumes from first principles with no free parameters, apart 
from the overall scale. 



GAUGE FIELDS ON THE THREE-SPHERE 

The reason to consider the three-sphere lies in the fact that the conformal equiv- 
alence of S 3 x 1R to H 4 allows one to construct instantons explicitly 11 ' 33 . This greatly 
simplifies the study of how to formulate 6 dependence in terms of boundary conditions 
on the fundamental domain, and indeed we will see that for S 3 simple enough results 
can be obtained to address this question 19 ' 34 . The disadvantage of the three-sphere is 
that in large volumes the corrections to the glueball masses are no longer exponential 9 . 

We will summarise the formalism that was developed in 11 . Alternative formula- 
tions, useful for diagonalising the Faddeev- Popov and fluctuation operators, were given 
in ref. 10 . We embed S 3 in IR 4 by considering the unit sphere parametrised by a unit 
vector It is particularly useful to introduce the unit quaternions er M and their 
conjugates a M = g^ by 

= {id, if), = (id, -if). (11) 

They satisfy the multiplication rules 

V^Gu = r]* v G a , G^G u =fj® u G a , (12) 

where we used the 't Hooft rj symbols 22 , generalised slightly to include a component 
symmetric in fi and v for a — 0. We can use r\ and fj to define orthonormal framings 35 
of S* 3 , which were motivated by the particularly simple form of the instanton vector 
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3.210 



> A 



Figure 5. "Artistic impression" of the potential for the three-torus. Shown are 
two vacuum valleys, related to each other by a gauge transformation h\ with wind- 
ing number 1, whose shape in the direc- tion perpendicular to the valley depends 
on the position along the valley. The induced one-loop effec- tive potential, of 
height 3.210/L, has degenerate minima related to each other by the anti-periodic 
gauge transformations g^)- The classical barrier, separating the two valleys, has a 
height 72.605/Lc/ 2 . 



potentials in these framings. The framing for S 3 is obtained from the framing of H 4 
by restricting in the following equation the four-index a to a three-index a (for a = 
one obtains the normal on S 3 ): 

^ = r\l v n v , el = fj^n v . (13) 

Note that e and e have opposite orientations. Each framing defines a differential oper- 
ator and associated (mutually commuting) angular momentum operators L\ and L^- 

L H S< • <"> 

It is easily seen that L\ = L 2 ,, which has eigenvalues 1(1 + 1), with I — 0, |, 1, ■ • •. 

The (anti-)instantons 36 in these framings, obtained from those on 1R 4 by inter- 
preting the radius in 1R 4 as the exponential of the time t in the geometry S 3 x 1R, 
become 

e-a ~r* o A e — (u + e ■ n)a . . 

A) = —, r, A = v-^-, (15) 

" 2(l + e-nY 2(1 + e-n) ' K ' 



where 



2s 2 2sb 

U = 1 i U2 \ 2 ' E V = 1 i 2 ' S = Ae< - ( 16 ) 

1 + b 2 + s 2 1 + b 2 + s 2 



Here e and A are defined with respect to the framing for instantons and with respect 
to the framing e° for anti-instantons. The instanton describes tunnelling from A = 
at t = — oo to A a = — a a at t = oo, over a potential barrier at t = that is lowest when 
6 M = 0. This configuration corresponds to a sphaleron 37 , i.e. the vector potential A a = 
— \a a is a saddle point of the energy functional with one unstable mode, corresponding 
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to the direction (u) of tunnelling. At t — oo, A a = —a a has zero energy and is a gauge 
copy of A a = by a gauge transformation g = n ■ a with winding number one. 

We will be concentrating our attention to the modes that are degenerate in energy 
to lowest order with the modes that describe tunnelling through the sphaleron and 
"anti-sphaleron" . The latter is a gauge copy by a gauge transformation g = n ■ a 
with winding number —1 of the sphaleron. The two dimensional space containing the 
tunnelling paths through these sphalerons is consequently parametrised by u and v 
through 

A^) = (-<-<)y- (17) 

The gauge transformation with winding number —1 is easily seen to map (u, v) = (w, 0) 
into (u,v) = (0,2 — w). The 18 dimensional space is defined by 

d ) = ( C K + d H) y = Mc > d)e > (18) 

The c and d modes are mutually orthogonal and satisfy the Coulomb gauge condition: 

d l A l (c,d) = 0. (19) 

This space contains the (u, v) plane through c" = — uS® and df = ~v5^. The significance 
of this 18 dimensional space is that the energy functional 11 

V(c, d) = - J s3 \ tr (i* ) = V(c) + V(d) + ^ {W) 2 (d*) a - (c^) 2 } , (20) 

V(c) = 2vr 2 |2(c«) 2 + 6 det c + |[(c?c?) 2 - (cfcj) 2 ]} , (21) 

is degenerate to second order in c and d. Indeed, the quadratic fluctuation operator At 
in the Coulomb gauge, defined by 

" L \ tr (jF; ^ = I S s tr (^«A0 + o (a 3 ) , 

At,, = 2L\5 l3 + 2 (Lx + Sf , ^. = -i£ aij , (22) 

has A(c, d) as its eigenspace for the (lowest) eigenvalue 4. These modes are consequently 
the equivalent of the zero-momentum modes on the torus, with the difference that their 
zero-point frequency does not vanish. 

FPf(A) in eq. (4) is defined as a hermitian operator acting on the vector space C 
of functions g over S 3 with values in the space of the quaternions 1H = {g M a M |g M G R}. 
The gauge group Q is contained in C by restricting to the unit quaternions: Q = {g G 
£\9 — 9^ a ^i9^ e R ) fl , Aifi , M — 1}- F° r arbitrary gauge groups C is defined as the algebra 
generated by the identity and the (anti-hermitian) generators of the algebra. When 
minimising the same functional over the larger space C one obviously should find a 
smaller space Ac A. Since £ is a linear space A can also be specified by the condition 
that FPf(A) be positive, 

A = {A G T| (g, FPf(A) g) > 0, V 5 G £}. (23) 

As for the Gribov horizon, the boundary of A is therefore determined by the location 
where the lowest eigenvalue vanishes. For the (c, d) space it can be shown 19 that 
the boundary dA will touch the Gribov horizon dQ. This establishes the existence 
of singular points on the boundary of the fundamental domain due to the inclusion 
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Figure 6. Location of the classical vacua (large dots), sphalerons (smaller dots), the Gri- 
bov horizon (fat sections), the boundary of A (dashed curves) and part of the boundary 
of the fundamental domain (full curves). Also indicated are the lines of equal potential 
in units of 2™ times the sphaleron energy. 

A C A C Q. By showing that the fourth order term in eq. (2) is positive (see app. A of 
ref. 19 ) this is seen to correspond to the situation as sketched in figure 2. 

To simplify the notation, write FPf = FPi and FP = FPi, with the indices 
related to the isospin. The associated generators are 

fx = if and T\ = iadf. (24) 

One can now make convenient use of the SU(2) 3 symmetry generated by L 1 , L 2 and T t 
to calculate explicitly the spectrum of FP t (A). One has 

FP t (A(c, d)) = AL\ - jC«T t a L\ - 2 -d«T?U 2 , (25) 

which commutes with L\ = L|, but for arbitrary (c, d) there are in general no other 
commuting operators (except for a charge conjugation symmetry when t = \). Re- 
stricting to the (u, v) plane one easily finds that 

FP t (A(u, v)) = AL\ + -uU ■ f t + -vL 2 ■ f t , (26) 

which also commutes with the total angular momentum J t = L\ + L 2 + T t and is easily 
diagonalized. Figure 6 summarises the results for this (u, v) plane and also shows the 
equal-potential lines as well as exhibiting the multiple vacua and the sphalerons. As it 
is easily seen that the two sphalerons are gauge copies (by a unit winding number gauge 
transformation) with equal norm, they lie on <9A, which can be extend by perturbing 
around these sphalerons 38 . 

To obtain the result for general (c, d) one can use the invariance under rotations 
generated by L\ and L 2 and under constant gauge transformations generated by T t , to 
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Figure 7. The fundamental domain (left) for constant gauge fields on S 3 , with respect to 
the "instanton" framing e" in the diagonal representation A a = x a a a (no sum over a). 
By the dots on the faces we indicate the sphalerons, whereas the dashed lines represent 
the symmetry axes of the tetrahedron. To the right we display the Gribov horizon, 
which encloses the fundamental domain, coinciding with it at the singular boundary 
points along the edges of the tetrahedron. 



bring c and d to a standard form, or express det ^FP ( (i(c, d))\ l= i j, which determines 

the locations of dQ and <9A, in terms of invariants. We define the matrices X and Y 
by X% = (cc*)£ and Y" b a = (dd 1 )^ which allows us to find 

det (FPi(A(c,d))\ l= i) = [81 - 18 Tr (X + Y) + 24(det c + det d) 

-(Tr(X-F)) 2 + 2Tr((X-F) 2 )] 2 . (27) 

The two-fold multiplicity is due to charge conjugation symmetry. The expression for 
t = 1, that determines the location of the Gribov horizon in the (c, d) space, is given 
in app. B of ref. 19 . If we restrict to d = the result simplifies considerably. In that 
case one can bring c to a diagonal form c" = XiSf. The rotational and gauge symmetry 
reduce to permutations of the Xi and simultaneous changes of the sign of two of the x,. 
One easily finds the invariant expression ( Tr (X) = J2i x \ an d det c = \[ i Xi) 

det (FPi(A(c,0))| J= i) = (2 det c - 3 Tr (X) + 27) 4 . (28) 

In figure 7 we present the results for A and Q. In this particular case, where d = 0, 
A coincides with A, a consequence of the convexity and the fact that both the sphalerons 
(indicated by the dots) and the edges of the tetrahedron lie on <9A, the latter also lying 
on dfl. It is essential that the sphalerons do not lie on the Gribov horizon and that 
the potential energy near dfl is relatively high. This is why we can take the boundary 
identifications near the sphalerons into account without having to worry about singular 
boundary points, as long as the energies of the low-lying states will be not much higher 
than the energy of the sphaleron. It allows one to study the glueball spectrum as a 
function of the CP violating angle 9, but more importantly it incorporates for 6 = 
the noticeable influence of the barrier crossings, i.e. of the instantons. 

An effective Hamiltonian for the c and d modes is derived from the one-loop effec- 
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Table 1. Numerical coefficients for V 



Kt = -0.2453459985179565 


K 2 = 


3.66869179814223 


k 3 = 0.500703203096610 


K4 = 


-0.839359633413003 


n 5 = -0.849965412245339 


K 6 = 


-0.06550330854836428 


n 7 = -0.3617122159967145 


K 8 = 


-2.295356861354712 



tive action 34 . To lowest order it is given by 

where g(R) is the running coupling constant (related to the MS running coupling by a 
finite renormalisation, such that kinetic term above has no corrections). The one- loop 
correction to the effective potential 34 is given by (for table 1): 

vS(c, d) = V<S(c) + V£(d) + «7(cf) 2 (dJ) 2 + ^ 8 (c^) 2 , (30) 
V$(c) = K^clf + k 2 det c + k 3 (cX? + « 4 (c?c^) 2 + /€ 5 (c^) 2 det c + k 6 (c^) 3 . 

Errors due to an adiabatic approximation are not necessarily suppressed by powers of 
the coupling constant. Nevertheless, one expects to achieve an approximate under- 
standing of the non-perturbative dynamics in this way 34 . 

The boundary conditions are chosen so as to coincide with the appropriate bound- 
ary conditions near the sphalerons, but such that the gauge and (left and right) rota- 
tional invariances are not destroyed. Projections on the irreducible representations of 
these symmetries turned out to be essential to reduce the size of the matrices to be 
diagonalised in a Rayleigh-Ritz analysis. Remarkably all this could be implemented 
in a tractable way 34 . Results are summarised in figure 8. One of the most important 
features is that the 0~ glueball is (slightly) lighter than the + in perturbation theory, 
but when including the effects of the boundary of the fundamental domain, setting 
in at / = g 2 /2n 2 ~ 0.2, the 0~/0 + mass ratio rapidly increases. Beyond / ~ 0.28 
it can be shown that the wave functionals start to feel parts of the boundary of the 
fundamental domain which the present calculation is not representing properly 34 . This 
value of / corresponds to a circumference of roughly 1.3 fm, when setting the scale as 
for the torus, assuming the scalar glueball mass in both geometries at this intermediate 
volume to coincide. 



DISCUSSION 

We have analysed in detail the boundary of the fundamental domain for SU(2) 
gauge theories on the three-torus and three-sphere. It is important to note that it 
is necessary to divide A by the set of all gauge transformations, including those that 
are homotopically non-trivial, to get the physical configuration space. All the non- 
trivial topology is then retrieved by the identifications of points on the boundary of the 
fundamental domain. 

As we already mentioned in the introduction, the knowledge of the boundary iden- 
tifications is important in the case that the wave functionals spread out in configuration 
space to such an extent that they become sensitive to these identifications. This hap- 
pens at large volumes, whereas at very small volumes the wave functional is localised 
around A = and one need not worry about these non-perturbative effects. That 
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0.1 0.25 

Figure 8. The full one-loop results for the masses of scalar, tensor and odd glueballs 
on S* 3 as a function of / = g 2 (R)/2ir 2 for 6 = 0. The dashed lines correspond to the 
perturbative result. 

these effects can be dramatic, even at relatively small volumes (above a tenth of a 
fermi across), was demonstrated for the case of the torus. However, for that case the 
structure of the fundamental domain (restricted to the abelian zero-energy modes) is a 
hypercube and deviates considerably from the fundamental domain of the three-sphere. 
Nevertheless, the spectrum for the sphere is compatible with that for the torus in vol- 
umes around one fermi across 39 , with m(2 + )/m(0 + ) ~ 1.5 and m(0~)/m(0 + ) ~ 1.7. 

It should be noted that the shape of A is independent of L if the gauge field is 
expressed in units of 1/L. Suppose that the coupling constant will grow without bound. 
This would make the potential irrelevant and makes the wave functional spread out over 
the whole of field space (which could be seen as a strong coupling expansion). If the 
kinetic term would have been trivial the wave functionals would be "plane waves" on 
a space with complicated boundary conditions. In that case it seems unavoidable that 
the infinite volume limit would depend on the geometry (like T 3 or S* 3 ) that is scaled-up 
to infinity. With the non-triviality of the kinetic term this conclusion cannot be readily 
made and our present understanding only allows comparison in volumes around one 
cubic fermi. However, one way to avoid this undesirable dependence on the geometry 
is that the vacuum is unstable against domain formation. As periodic subdivisions 
are space filling on a torus, this seems to be the preferred geometry to study domain 
formation. In a naive way it will give the correct string tension (flux conservation tells 
us to "string" the domains that carry electric flux) and tensor to scalar mass ratio 
(averaging over the orientations of the domains is expected to lead to a multiplicity 
weighted average of the T 2 and E masses). Furthermore, the natural dislocations of 
such a domain picture are gauge dislocations. The point-like gauge dislocations in four 
dimensions are instantons and in three dimensions they are monopoles. Their density 
is expected to be given roughly as one per domain (with a volume of around 0.5 cubic 
fermi). Also the coupling constant will stop running at the scale of the domain size. 
We have discussed this elsewhere and refer the reader to refs. 6 ' 32 ' 40 for further details. 
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